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While much of the focus of observational cosmology is on the high redshift universe it is important not to 
neglect the very local universe as a source of cosmological information. The inner profiles and number counts 
of low mass halos have provided the biggest stumbling blocks so far for ACDM. These small structures can 
only be seen nearby In the very local universe {cz < 3000 km s~^) the component of a galaxy’s redshift 
which is due to motions under gravity can be comparable to (or even larger than) its cosmological redshift. 
The distance to a galaxy as inferred from its redshift can differ by more than a factor of two from its actual 
distance. Given that the mass and intrinsic size scale of a galaxy (among many other physical parameters) 
depend strongly on distance, and that these peculiar motions are coherent over large regions of the sky, serious 
biases can occur. While it is important to have accurate distances to local galaxies, it is not feasible to measure 
a primary distance for every nearby galaxy. Instead, a velocity field model can be used to provide a first order 
correction to the redshift distance. We report on a new sample of Tully-Fisher distances (the SFlH—h) which 
is being used in combination with publicly available primary distances to model galaxy flows in and around 
the Local Supercluster. This sample has '^lO times as many tracers as were used for the current best model. 
Initially a parametric model including infall onto multiple attractors will be used. Such models (which assume 
spherical symmetry for the attractors) are not realistic representations of the true velocity field, but provide 
useful first order corrections. We explore the typical errors in distances from such a flow model and how they 
vary in different regions of the local supercluster by fitting the same parametric model to mock catalogs derived 
from a constrained simulation of the volume. Non-parametric reconstructions of the velocity field will follow. 


1. Introduction 

Cosmology has entered a new era. Recently large 
redshift surveys have been used in combination with 
observations of the cosmic microwave background to 
provide such strong constraints on the cosmological 
parameters that many astronomers now consider cos¬ 
mology ‘solved’ (eg. 0|), giving the concordance 
model of ACDM in which the universe consists of 
4% baryons, 26% dark matter and 70% dark energy. 
Much of this information has come from studies of the 
very high redshift universe - for example the WMAP 
observations of the fluctuations in the cosmic mi¬ 
crowave background Q, or the studies of high red¬ 
shift Type la supernova that gave the first evidence 
for dark energy (lllll . 0 

On the large scales ACDM has stood up extremely 
well to observational tests. Its biggest challenges have 
so far come from much smaller scales; for example in 
measurements of the inner density profiles and num¬ 
ber counts of low mass galaxies. Such observations 
can only be done in the local universe, and therefore 
are impacted by the uncertainty inherent in estimat¬ 
ing distances in this region where Hubble’s law cannot 
be used. There are many thousands of galaxies with 
measured redshifts less than a few thousand km s“^, 
so it is obviously not feasible to measure redshift in¬ 
dependent distances to all of them. Flow models are 
therefore the only was to provide unbiased distances 
to these galaxies. 

In this paper we will discuss the impact the local pe¬ 
culiar velocity field has on cosmology, and the use of 
simple flow models to estimate galaxy distances in the 
local universe. In Section 2 we will describe the local 
velocity field and the pattern of deviations from Hub¬ 


ble’s Law. In Section 3 we will discuss an example of 
the impact of ignoring peculiar velocities, specifically 
on the inferred number counts of low mass galaxies. In 
Section 4 we will introduce the SFlH—h, a new sample 
designed to study the local velocity field. In Section 5 
we provide a preliminary report on a test of distances 
estimated using simple infall models, which is done 
by fitting such a model to a mock catalog generated 
from a constrained simulation of the volume. Section 
6 provides a summary of the paper. 


2. The Local Velocity Field 


The uneven distribution of mass in the universe 
gives galaxies deviations from Hubble flow (or peculiar 
velocities) of the order 200-600 km s“^. In the local 
universe cosmological redshifts are small, so these pe¬ 
culiar velocities can dominate the observed recessional 
velocity of a galaxy. Gravitational instability theory 
in the linear regime in an expanding universe shows 
that the peculiar velocity at any given point in space 
is directly proportional to the gravitational potential. 


v(x) 


2 


/g(x) 


( 1 ) 


where the proportionality constant depends on cos¬ 
mological parameters, and the factor / ~ which 
describes the rate of growth of structures. N-body 
simulations have shown that these motions are domi¬ 
nated by the infall of galaxies along filaments towards 
the clusters which form at the intersections of fila¬ 
ments. 

Our nearest large cluster is the Virgo Cluster at a 
distance of ~16 Mpc away, and within the Local Su- 
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Figure 1: A slice of the density field in the Supergalactic plane from the constrained simulations described in [13 . The 
density field has been smoothed with a Gaussian of radius 4 Mpc. Contours are spaced at intervals of 
5 = (p — po)/po = 0.2. The mean density of the universe is shown by the thick solid line, while the magenta solid lines 
show overdensities, blue dotted lines underdensities. 


percluster (of size ^30 Mpc) galaxies lie preferentially 
in a plane, known as the supergalactic planeQ. At 
our location the infall velocity onto Virgo is Uinf ~200 
km s“^, but this infall cannot account for all of the 
V = 368±1.8 km s“^ peculiar velocity of the Local 
Group in the CMB frame (as measured by COBE 
and now WMAP Q), largely because it points in a 
direction '^50° away from the direction to the Virgo 
cluster. This and other observations of the motions of 
galaxies in the local universe have revealed a flow in 
a the direction towards the Hydra-Centaurus cluster 
and Shapley Supercluster in the background (or the 
“Great Attractor” region), as first suggested by Eg. 

The distribution of mass in a slice through the su¬ 
pergalactic plane from the constrained simulations de¬ 
scribed in E2 shown in Fig ^ with the major 
nearby concentrations of mass indicated. This slice 
goes out in distance to cz = 8000 km s“^ and has been 
smoothed with a Gaussian of radius 4 Mpc {Hq = 70 
km s“^ Mpc“^ in the simulation). 

3. The Missing Sateiiite Probiem 

One of the largest discrepancies between observa¬ 
tions and AGDM revolves around the “missing satel¬ 
lite” problem. Numerical simulations in ACDM pre¬ 


dict a value for the logarithmic slope of the faint end 
of the halo mass function, a ~ —1.8, while most recent 
determinations of the optical luminosity function (LF) 
in the local volume yield values of a that are signif¬ 
icantly flatter. Gomplementary to the results on the 
optical LF, several determinations of the local HI mass 
function (HIMF) have likewise yielded relatively flat 
faint-end slopes. The over-prediction of the number of 
low mass objects relative to those actually observed is 
considered one of the last remaining stumbling blocks 
for ACDM, and points towards a new understanding 
of the baryon physics of galaxy formation. It is ob¬ 
viously desirable that all of the uncertainties in the 
observed luminosity and mass functions be well un¬ 
derstood. 

To derive a luminosity (or mass) function from ob¬ 
servational data, the sample selection effects must be 
considered carefully. In most surveys, the lowest lu¬ 
minosity (or mass) objects are only visible nearby, 
and their contribution must be weighted accordingly. 
Fractional distance errors can be large for relatively 
nearby objects, and can thus have a very strong im¬ 
pact on determinations of the faint end of luminosity 
(or mass) functions. In particular, the derivations of 
HIMFs have to rely on relatively shallow, wide area 
surveys to detect low HI mass objects over a sufficient 
volume, and are therefore most susceptible to uncer- 
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Figure 2: Difference between input and reconstructed HIMFs for a mock HIPASS BGC survey (shown are the 
average of 10 simulations). The HIMF derived when “true” distances are used (blue circular points and dotted line), 
along with the HIMF constructed when multi-attractor model distances are used (green triangular points and 
long-dashed line) reconstructs the input HIMF within the errors (here typical Poisson counting errors are shown only 
on the “true” distance points). When pure Hubble flow is assumed (red square points and short dashed line), the low 
mass end of the HIMF is underestimated. The solid line shows the HIMF derived in 1^ . which differs from the input 
by Aa = 0.1. 


tainties associated with distance errors in the nearby 
universe. 

In El we described simulations of the impact of ne¬ 
glecting peculiar velocities on the derivation of HIMFs 
from recent HI surveys. We showed that the HIMF 
derived from the HIPASS Brightest Galaxy Catalog 
underestimated the slope of the low mass end of 
the HIMF because of their use of Hubble’s law to esti¬ 
mate distances to nearby galaxies (see Figl^J. HIPASS 
is a HI survey of the southern sky. Nearby galaxies in 
the south will systematically have their distances (and 
therefore masses for a given HI flux) overestimated 
if peculiar velocities are neglected, largely because of 
their motion towards the “Great Attractor” region at 
(RA, DEC) ^ (13 hr, -35°). This systematic effect 
results in an underestimate of the number counts of 
low mass galaxies as inferred from the HIMF derived 
from this data. 

The size of the bias is not nearly large enough to 
explain the missing satellite problem, but it serves to 
illustrate the point that without a good understand¬ 
ing of the local peculiar velocity field mistakes can be 
made which influence our understanding of cosmology 
as a whole. 


4. The SFI++ 

The SFI-I--I- is a sample of spiral galaxies with I- 
band photometry and a combination of HI and Ha 
spectroscopy suitable for use in the Tully-Fisher (TF) 
relation (a relation between the absolute magnitude 


and rotation velocity of a spiral galaxy). This sample 
builds on the all-sky SFI (Spiral Field Tband) sample 
(0)01)0 10 ) contains much new I-band photom¬ 
etry, HI and Ha spectroscopy in areas of the sky north 
of -20°. Data from the SFI has also been reprocessed 
to provide a uniform sample. The SFI-I--I- will pro¬ 
vide Tully-Fisher distances for ~ 5000 galaxies out to 
cz = 10 000 km s“^. It is a diameter limited catalog, 
with diameter limits which vary with the redshift of 
the galaxy in a similar way to the SFI, except that the 
minimum diameter limit for SFI-I—h is smaller. Fig¬ 
ure 13 shows the preliminary sky distribution of the 
galaxies in this sample. 

A multiattractor model is being fit to a subset 
of these Tully-Fisher distances (for galaxies within 
cz =4000 km s“^) combined with a sample of ^ 500 
publicly available primary distances collected from the 
literature. Primary distances provide lower errors for 
individual measurements, but larger numbers of ve¬ 
locity held tracers can be obtained using secondary 
indicators like TF. This combination will provide the 
best currently available sample to study the local ve¬ 
locity held and represents a signihcant increase in the 
number of velocity held tracers from what has pre¬ 
viously been used. The multiattractor model should 
not to be considered a realistic representation of the 
local density and velocity helds, but will provide hrst 
order corrections to galaxy distances estimated from 
their redshifts. Non-parametric modeling of the ve¬ 
locity held (e.g. POTENT j^) which can provide a 
more realistic picture will follow. 
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Right Ascension 


Figure 3: The sky distribution of galaxies in the SFI++ sample centered at RA=12 hours. Galaxies are shown as blue 
filled circles if they have a negative peculiar velocity, red open circles if they have a positive peculiar velocity. The lines 
show the location of the Galactic plane. 


5. Testing Multiattractor Models on 
Constrained Simulations 

Numerical simulations in the ACDM paradigm and 
redshift surveys have given us a consistent picture of 
the filamentary large scale structure of the universe 
in which the dominant motions of galaxies are that 
of infall onto clusters along filaments. In light of this 
picture the traditional infall models used to study the 
velocity field of the Local Supercluster going back to 
[l^. and most recently used by ^.t 

best over simplifications and at worst simply wrong. 
The assumption of spherical symmetry for attractors 
is clearly an over simplification, even for a single clus¬ 
ter like Virgo which is made up of many sub-groups, 
let alone for a filamentary or pancake like superclus¬ 
ter. 

However, multiattractor models have many nice fea¬ 
tures which leads to their continuing use. They are 
simple to apply, and the non-linear approximation 
of 1^ (which is only valid for spherical masses) has 
been shown to work reasonably well to overdensites 
of (5 ^ 30 0. For a galaxy which has only a posi¬ 
tion on the sky and a redshift, a multiattractor model 
will provide a unique distance over much (if not all) 
of its volume. A multiattractor model can in principle 
have any number of spherical attractors which could 
be combined to give more complicated mass distri¬ 
butions, and other components, such as quadrupole 
moments in the velocity field (perhaps from the tidal 
influence of distant masses) are also easy to add. A 


multiattractor model can also give you an idea of the 
dominant sources of attraction in the local volume, 
taking into account a combination of the total mass 
and distance to the attractors. 

Here we test the validity of a simple multiattractor 
model as applied to a mock distance catalog gener¬ 
ated from the Virgo Consortium’s publicly available 
constrained simulations of the local universe da We 
test distance predicted from the best fit multiattrac¬ 
tor model against the real distances in the constrained 
simulation to see how well such a model can be used 
to predict distances for galaxies in the local universe. 


5.1. Mock Catalog of Distances 

The mock catalogs used here are based on the con¬ 
strained simulations described in na These ACDM 
simulations are designed to mimic the local densiW 
and velocity field. Semi-analytic galaxy formation Q 
is grafted onto the N-body simulation to create galaxy 
catalogs. We generate simple mock catalogs of dis¬ 
tance indicators from these publicly available catalogs 
by imposing a magnitude limit of mj = 12 and se¬ 
lecting a random sample of ~ 500 such galaxies with 
cz < 4000km s“^. The sample is uniform over the 
sky except that no galaxies are seen within 10° of the 
Galactic plane. Gaussian scatter of width 0.2 mag is 
added to the distance modulus to simulate 10% dis¬ 
tance errors. 
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Figure 4: Contours of constant heliocentric velocity in the snpergalactic plane for the best fit multiattractor model of a 
mock sample of distance indicators. The location of the mock galaxies used in the fit are also shown. Thin black 
contours are places at 200 km s“^ intervals with the thick red contours every 2000 km s“^. The concentric green 
dashed circles show pure Hubble flow with Ho = 70 km s“^ Mpc“^. 


5.2. The Best Fit Multiattractor Model 

We fit a multiattractor model to the mock distance 
catalog described above. This model has a very simi¬ 
lar functional form to that described in 0. except 
that we use NFW profiles for the attractors. The 
best fit model includes infall onto both Virgo, and 
the “Great Attractor” region and a quadrupole com¬ 
ponent to the velocity field which might account for 
the tidal influence of external masses. Fig 0 shows 
contours of constant heliocentric velocity in the su- 
pergalactic plane from this model, along with the lo¬ 
cations of the mock galaxies it was fit to. 


5.3. Predicting Distances from a Flow 
Model 

A flow model gives you the peculiar velocity of a 
galaxy at a given location. In order to use it to give 
distances for galaxies with observed redshifts one must 
invert it, and therefore solve an equation of the form 
D = /(£>, nobs, R-A,DEC). We do this using Numer¬ 
ical Recipes root finding algorithms. Galaxies which 
are close to the location of attractors are assumed to 
be in the attractor core. Galaxies in the triple values 
regions (where three distances along a line of sight 
have the same observed recessional velocity) are as¬ 
signed to the central location, and the distance error 
is inflated to bracket all three possible distances. 


5.4. Flow Model Distance Errors 

The use of a this flow model gives distance esti¬ 
mates to galaxies to less than 5% accuracy over much 
of the local volume. This is illustrated in Fig. 0which 
shows contours of the average percent error for dis¬ 
tances estimated using the flow model. Fig. |H1 shows 
the same thing, but this time for distances estimated 
using Hubble’s Law. Note that near the origin where 
the absolute distance is small there will always be an 
area where distances are overestimated by large per¬ 
centages (it is not possible to underestimate distances 
to galaxies which are very close by). The flow model 
does not completeley correct for the area of systemat¬ 
ically overpredicted distances in the direction towards 
the “Great Attractor” region at (SGX, SGY) ~ (-50, 
20) Mpc, nor can it completely correct for the velocity 
field in the direction towards Virgo at (SGX, SGY) ~ 
(-3, 20) Mpc. 

6. Summary 

We have argued that the local universe is an impor¬ 
tant source for cosmological information, especially in 
light of the fact that the two largest challenges still 
faced by AGDM come from observations which are 
only possible in the local universe. 

Within a redshift of cz ~ 4000km s“^ Hubble’s law 
cannot be used to predict distances. We show, using 
the example of the derivation of the HI mass function 
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Figure 5: Contours of the mean percent error in 
distances estimated using the multiattractor model for 
galaxies within 5 Mpc of the supergalactic plane. Errors 
are smoothed on a scale of 4 Mpc. The contour spacing 
shows a 5% error with red solid contours being 
overestimated distances, blue dahsed contours 
underestimated. 



Figure 6: As in Fig |21 except for distances estimated 
assuming Hubble’s Law. 

for galaxies in the southern sky, that if Hubble’s law is 
used to derive distance serious biases can occur. This 
work is presented in detail in 0. 

We present the SFI++, and all-sky sample of ^5000 
spiral galaxies for which Tully-Fisher distances will be 
available. This sample builds on the SFI sample of the 
1990s (0101010) and provides the best currently 
available sample for studies of the local peculiar ve¬ 
locity field. Initially a multiattractor model will be fit 
to this data, more realistic modeling will follow. 

We provide a preliminary report of work testing 
the distances predicted the using commonly applied, 
but simplistic multiattractor infall models. This work 


makes use of the publicl y av ailable Virgo Consortium 
constrained simulations |l^ , and fits a multiattractor 
model to a mock sample of 500 galaxies with distances 
measured to 10% derived from this simulation. The 
best fit model includes infall onto both the Virgo clus¬ 
ter, and Hydra-Centaurus (‘Great Attractor’) region, 
as well as a small quadrupolar component to the ve¬ 
locity field. Preliminary results suggest that this mul¬ 
tiattractor model can be used to estimate distances to 
better than 5% on average over much of the local vol¬ 
ume, but care must be taken in the directions towards 
both attractors. 
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